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CHAPTER 1 


INTRODUCTION 


1.1 OZONE MEASUREMENT 

The measurement of ozone in the atmosphere has become increasingly important 
over the past two decades. Significant increases of ozone concentrations in the lower 
atmosphere, or troposhpere, and decreases in the upper atmosphere, or stratosphere, have 
been attributed to man-made causes [1]. High ozone concentrations in the troposphere 
pose a health hazard to plants and animals and can add to global warming [1], On the 
other hand, ozone in the stratosphere serves as a protective barrier against strong 
ultraviolet (UV) radiation from the sun [2], Man-made CFC's (chloroflourocarbons) act 
as a catalyst with a free oxygen atom and an ozone molecule to produce two oxygen 
molecules therefore depleting the protective layer of ozone in the stratosphere [3], The 
beneficial and harmful effects of ozone require the study of ozone creation and 
destruction processes in the atmosphere. Therefore, to provide an accurate model of 
these processes, an ozone lidar system must be able to be used frequently with as large a 
measurement range as possible [4], 

Various methods can be used to measure atmospheric ozone concentrations. 
These include different airborne and balloon in situ measurements, solar occulation 
satellite techniques, and the use of lasers in lidar (light detection and ranging) systems to 
probe the atmosphere [4], Typical in situ devices such as weather balloons can only 
measure within the direct vicinity of the instrument and are therefore used infrequently 
[4], Satellites use solar occulation techniques that yield low horizontal and vertical 
resolution column densities of ozone. Lidar uses lasers to probe the atmosphere and 



analyze the backscattered laser energy to yield information such as temperature 
distributions and molecular densities. Lasers provide a means in which long-range non- 
intrusive measurements of ozone as a function of altitude can be made with very high 
vertical resolution. The range of lidar systems extends to tens of kilometers making 
ground-based measurements of the high altitude stratosphere possible. 

1.2 LIDAR 

Lidar has been in use since 1963 [5,6] and works on the same principal as 
ordinary radar. A radar uses a transmitter to send a directed pulse of electromagnetic 
radiation into a medium such as air. If the pulse encounters an object in the medium then 
some of the energy is scattered or reflected. A radar receiver collects some of this 
reflected radiation, which contains information about the object that it has encountered. 
For example, the delay between the transmission and detection of the reflected signal can 
be transformed into the distance from the radar transmitter to the object. Also, a moving 
object causes a frequency shift in the radar pulse. The speed of the object can be 
calculated from this shift. Lidar is most commonly used for remote sensing of the 
atmosphere. One common use of lidar is to measure the concentration of different types 
of molecular species that are present in the atmosphere. 

A lidar system for atmospheric studies consists of a pulsed laser source, a 
telescope receiver to collect returned light, a light detector, and computer for data 
analysis. Lasers offer the advantage of coherent and collimated sources that allow long 
range remote sensing with high temporal and spatial resolution [5]. A laser is pulsed 
into the atmosphere interacting with the molecules and atoms present. A receiver 



telescope then collects the backscattered radiation, which in typical configurations is 
close to the laser source. The receiver telescope focuses the light return onto an optical 
detector, such as an avalanche photodiode (APD) or a photomultiplier tube (PMT), turns 
the returned light pulse into an electrical pulse. This detector is connected to other 
electronics such as amplifiers and filters and then sent to a computer for analysis. 

1.3 OZONE Differential Absorption Lidar (DIAL) 

Lidar cannot be used to measure ozone concentrations as a function of altitude, 
thus we use a technique called Differential Absorption Lidar (DIAL) to directly measure 
ozone concentrations. In a typical DIAL system, two laser pulses, separated in time and 
wavelength, are emitted by the laser system. These pulses, one tuned to high ozone 
absorption (on-line) and the other tuned to a lower absorption (off-line), are transmitted 
into and then backscattered by the atmosphere creating a return light signal. This signal 
is received by a telescope, passed through a narrow band filter, and focused onto a 
photomultiplier detector. When ozone is present in the atmosphere, the on-line signal 
decay is faster than that of the off-line due to the absorption by ozone in the atmosphere. 
The main advantage of the DIAL technique is that the number density of a species such 
as ozone can be deduced as a function of altitude. If the two lasers are close in 
wavelength then each laser pulse should undergo the same aerosol and molecular 
scattering. The main difference between the two pulses should only be the amount of 
absorption by ozone. If the absorption cross-sections are well known at each wavelength 
then accurate ozone number densities as a function of altitude can be calculated. 



1.4 OZONE DIAL SYSTEM AT NASA LANGLEY RESEARCH CENTER 

The Lidar Applications Group at NASA Langley Research Center has been 
successfully performing airborne DIAL measurements of tropospheric and stratospheric 
ozone and aerosols since 1980 [7], Many of the missions over the past eighteen years 
have been international experiments including missions over rain forests in Brazil during 
the burning season and missions to the Arctic to study the evolution of the ozone hole 
[Richter], The DIAL system has been continually upgraded over the years improving 
laser quality, receiver optics, detection electronics, and data analysis to increase the 
quality, accuracy, and range of atmospheric ozone and aerosol measurements [7], 

1.5 RESEARCH GOAL 

The goal of this research is to construct and test a complete lidar receiver system 
to measure stratospheric ozone concentrations. The system will include a fiber-optic 
coupled lidar telescope to collect the backscattered lidar signal and focus it onto a 
photomultiplier tube (PMT). Photon counting will be used to detect the weak signal 
return from the high altitude stratosphere. This receiver system will eventually be 
incorporated into the aircraft based instrument. The characteristics of the telescope 
including total transmission efficiency and limits of its measurement range will be 
determined. Also, optimum conditions for use of this telescope with a photon counting 
system will be made. The lidar equation will be solved for expected lidar return signals 
and these predictions will be compared to actual lidar returns. Finally, stratospheric 
ozone concentrations will be measured and compared to ozone sonde balloon 


measurements. 



CHAPTER 2 


THEORY 

This chapter covers the theory, derivation of equations, and a light detection 
technique used in lidar systems. First, the chapter begins with the derivation of the 
elastic backscattered lidar equation with emphasis on the effects of the geometry of the 
receiver telescope. The DIAL equation is then derived to determine the number density 
of an absorbing atmospheric species. Finally, photon counting and its advantages for 
lidar signal detection is covered. 

2.1 ELASTIC BACKSCATTERED LIDAR EQUATION 

The scattering form of the lidar equation yields the expected received power from 
a laser propagated into the atmosphere and collected by a receiver that is located near the 
laser system. The intensity variation in a lidar return signal can be used to determine the 
relative distribution of molecules and aerosols in the atmosphere. When a high density of 
molecules or aerosols is present, then the laser light is strongly backscattered and 
attenuated rapidly. This results in a strong initial signal with a fast decay rate. 
Conversely, when a low number of particles are present, the laser light propagates further 
into the atmosphere due to the smaller attenuation and scattering effects. The 
backscattered power received by the system P r (A,R) is 
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where P 0 is the initial power of the laser pulse, A (A,./?) is the receiver system function, R 
is the range from the receiver to the probed volume, c is the speed of light, r L is the laser 
pulse duration, /?(A,i?)is the cross section for scattering in the backward direction, and 

k(A,R) is the total extinction coefficient that accounts for all scattering and absorption 
losses. [4,5,6, 8], 

The receiver system function A(A,R) can become quite complex depending on 
the geometry of the optics in the receiver system. The system function must account for 
the transmission of the receiver optics, the spatial distribution of laser energy, and the 
range dependent overlap of the telescope's FOV and the laser beam (see Figure 2.1) [6]. 


d = d 0 -R8 



Figure 2.1 Overlap of telescope field of view and laser beam divergence. 
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Most modem lidar systems are monostatic in that the receiver telescope is located near 
the laser transmitter [6], In general there are two monostatic arrangements known as 
coaxial and biaxial. A coaxial arrangement has the telescope field of view axis 
coincident with the axis of the laser beam. This arrangement allows full overlap of the 
laser beam and the telescope's field of view at all ranges. In a biaxial system, the laser 
beam enters the field of view at some predetermined range. This arrangement allows 
discrimination against strong near-field signal returns that may saturate the light detector 
[4,5,6,7] (see Appendix A). 

The system function for a receiver separated by a distant d 0 from the laser 
transmitter then becomes 


A(A,R) = - £ ' ° j \^R,r,y/) F(R,r,y/)-r dr dy/ (2.2) 

\ iX J r - 0 y /~ 0 

where c{/Q accounts for the transmission of the optics and detector efficiency at 
wavelength a , A 0 is the area of telescope's primary mirror, W(R) is the radius the laser 
beam at altitude R , is the probability of backscattered light from point 

( r >^) (polar coordinates) being collected by the receiver, ./?,/•, ^)is the laser energy 
distribution at point (r,y/). and rj is the radius of the telescope's field of view. For a 
Gaussian laser beam the energy distribution FiyR,r,y/) is 



F{R, r ,y / ) = exp 
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In a monostatic biaxial system, the telescope FOV partially overlaps the laser beam at 
some altitudes. The partial overlap of the FOV and the laser at altitude R can be 

calculated as the intersection of two circles [6], The probability <^R,r,y/) of the 

backscattered light being collected can then be described as the area A (see fig.2. 1) of this 
intersection [6], This overlap function is range dependent with three possible situations 
occurring. 

1. At low altitudes the circles do not overlap, therefore A = 0. 

2. At high altitudes, the circles completely overlap and A is equal to the area of the 
smaller circle, which is usually the laser beam so A= n- W 2 (/?). 

3. At middle latitudes, the circles partially overlap. The area A is then equal to the area 
of the intersection of the two circles [6], 

A (r T ,W,d) = W 1 x w + r ]Xr -r T d s\n{x r ) (2.4a) 

where 


J* = cos 


+W^ -T£ 
2 ■ r T d 


(2.4b) 


Jr = COS 


c/~ +r[ -W^ 
2 • r T d 


(2.4c) 


Substituting Equations (2.3) and (2.4) into (2.2) may seem to make the integral term too 
complex, but this integral can be solved quite easily with computers using numerical 
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methods. This will be seen in Chapter 5 "Results" where the predicted signals, using the 
preceding theory, and actual signal returns are compared. 

2.2 DIFFERENTIAL ABSORPTION LEDAR (DIAL) 

Although the relative distribution of molecules and aerosols can be determined 
from a single lidar return, it is impossible to determine the absolute range resolved 
number density of a species in the probed volume. To determine the number density, the 
DIAL technique can be used to isolate the absorption effects of a particular molecular 
species. The DIAL method offers a very sensitive method to remotely measure the 
number density as a function of altitude. 

The DIAL method measures two individual lidar returns that have slightly 
different wavelengths. One of the wavelengths is tuned to a high absorption line (online) 
of the species and the other is tuned where less absorption occurs (offline). Since the 
wavelengths are only slightly different, each signal undergoes the same amount of 
molecular and aerosol scattering and absorption. The only difference in the return signals 
should come from the absorption by the species of interest, ozone in the present study. 
Figure 2.2 shows the absorption cross-section of ozone as a function of wavelength. The 
online and offline wavelengths transmitted by the NASA UV DIAL system are 300 and 
311 nm respectively. This results in absorption cross-sections from Figure 2.2 to be 
30x1 O ' 20 and 9x1 O ' 20 cm 2 / molec. 
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Wavelength (nm) 


Figure 2.2 Ozone absorption cross-section vs. wavelength. Cross-section values at 

wavelengths of 300 and 311 nm. 

The general form of the DIAL equation can be derived by taking the ratio of the 
return signals. The ratio of the online and offline return signals over a finite range 
between Rj and R 2 yields 


4 u .*) p.(k„r) 
" A(V«) 


exp 
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The system functions A (/t,/?) and the J R 2 dependence from each signal cancel out in the 
ratio leaving the backscatter and total extinction terms [4,6], The absorption feature of 




the species of interest such as ozone can be separated from the total extinction coefficient 


as 


ic(A t R) = jc a (A f R) + <7(A)Ar(R) (2.6) 

where K a is the absoprtion and extinction due to aerosols, a(A) is the ozone absorption 
cross section and N(R) is the ozone number density. If the absorption cross sections at 
the two wavelengths are well known, then we can solve for N(R). Defining the difference 
in absorption cross sections of the online and offline wavelengths as Act = <j on - <j off , 

gives the general form of the DIAL equation over a range cell A R = (R 2 -R,). 




1 

~ 2 -Aa AR ; 
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The [M] term yields the more familiar DIAL equation when the [B] and [E] terms are 
negligible [4,6], The [B] and [E] terms account for differences in the backscatter and 
extinction that result from the two different on and off wavelengths. If the difference in 
wavelength is small, the backscatter is spatially homogeneous, and the sky is relatively 
clear (high visibility, no clouds), then both terms are small and can be ignored [Browell]. 
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2.3 PHOTON COUNTING SYSTEMS 

Photon counting is a common light detection technique used in lidar applications. 
This section begins with the basic operating characteristics of photomultiplier tubes as 
light detectors used in photon counting systems. The theory of photon counting and the 
components of a typical system will be discussed. Also, the different sources of noise 
and precautions to reduce noise will be covered. Finally, a brief derivation of the signal- 
to-noise ratio will follow. 

2.3.1 Photomultiplier Tubes(PMT) 

Photomultiplier tubes (PMTs) are very sensitive light detectors. A typical PMT 
has three major components: a photocathode, a dynode chain, and an anode (see Fig. 2.3). 
The photocathode is a photosensitive material coated on the internal surface of the PMT 
that releases electrons when light strikes it. A high negative voltage is placed across the 
photocathode and dynode chain through a series of resistors. The electrons released from 
the photocathode are accelerated and focused by the high electric field toward a CsSb 
metal target known as the dynode. The dynode releases many electrons when one 
electron strikes the surface. These secondary electrons pass through a chain of dynodes 
amplifying the number of electrons as they hit each dynode. At the end of the chain, the 
electrical pulse collected by the anode can contain as many as 10 8 electrons [9], 




Figure 2.3, Typical photomultiplier tube 

Two important characteristics of PMTs are the quantum efficiency (QE) and the gain. 
The QE describes the number of electrons released for each incident photon. For 
example, if the QE is 20%, then 20 electrons should be released for every 100 photons 
incident on the PMT. It is desirable to have as high a QE as possible in order to have a 
high signal-to-noise ratio. The gain of a PMT is the total number of electrons that arrive 
at the anode from one electron released from the photocathode. The electron is focused 
on the first dynode that has a secondary electron emission coefficient of 5. In a chain of 
n dynodes, the total number of electrons at the anode will be 5 n . 


2.3.2 Photon Counting theory 

When a photon falls on the photosensitive cathode of the PMT, an electron is 
released and accelerated toward the first dynode of the dynode chain. The chain 
increases the number of electrons and a pulse of electrons is collected at the anode. In 
normal analog detection mode, the anode is connected to a load resistor and the electron 
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pulses create an average DC current. Typically, the voltage signal across the load resistor 
is connected to an amplifier and then digitized and averaged for computer analysis. The 
analog detection method serves as a good detection technique when high light intensities 
are being measured because of the high signal-to-noise ratio. But, when the light 
intensity is low, the noise in the system can be higher than the signal putting a limit on 
the sensitivity of the PMT detector. For DIAL applications, the lowest detectable signal 
translates into the maximum atmospheric measurement range of the system. 

For low light intensity measurements, the discrete nature of the photon can be 
used. The output pulses at the anode of a PMT are fast (several nanoseconds), separated 
in time, and discrete. Hence, these fast pulses can be digitally counted instead of 
measured as an average current. This type of PMT operation is known as single electron 
response (SER), which describes the PMT's ability to create an electron pulse from a 
photoelectron released from the photocathode [9], Photon counting with this type of 
PMT operation can offer significant improvements in signal-to-noise ratio compared to 
analog measurement techniques. 

A typical photon counting system, shown in Figure 2.4, consists of a PMT, an 
amplifier, a pulse height discriminator, and a digital counter. Photons from the light 
source strike the photocathode producing electrical pulses at the anode. The output 
pulses from the PMT are sent through a high bandwidth amplifier. The amplifier brings 
the PMT signal to a level compatible with the pulse height discriminator. The pulse 
height discriminator compares the input signal to an adjustable threshold voltage. All 
input pulses above this threshold are sent to the output with equal pulse height and width 
so they can be digitally counted. Any pulses below the set threshold are ignored and do 
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not get sent to the output. The digital counter is a multi-channel scaler with an averaging 
memory. Each channel, or bin, stores the number of pulses received from the 
discriminator during a count time determined by the system clock. For example, for a 
common clock rate of 1 MHz, each bin stores the number of pulses received in 1 usee. 




Figure 2.4 A typical photon counting system 



The output pulse of electrons at the anode is characterized by its pulse height. 
The pulse height is determined by the total number of electrons that reach the anode (see 
Figure 2.5). 


O ► 
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SINGLE 
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O ► 


ELECTRON 
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Figure 2.5 Pulse height of electrical output pulse of a PMT 

Variations in pulse height result from uncertainties in the secondary electron emission 
coefficient of the dynodes and the spread in electron trajectories and velocities [9], Some 
of the electrons deviate from their normal trajectories and therefore do not participate in 
the multiplication process. Other variations in pulse height are caused by noise pulses 
that do not originate from the photocathode. For example, electrons may be emitted from 
the dynodes by thermionic emission. The pulse height that results from these electrons is 
smaller because they do not travel through the entire multiplication chain. 

The pulse height distribution (PHD) of a typical PMT with and without a light 
signal applied is shown in Figure 2.6. This PHD shows that the bulk of the noise, which 
is mostly from dark current, is in the lower pulse height region. The PHD of the applied 
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light signal plus the noise reveals that pulses from the signal to be measured are mainly 
distributed in a finite region known as the single electron response region. 

The optimum measurement point of the PHD is the single electron response peak 
shown in Figure 2.6. The discriminator is set at this peak to ignore pulses below this 
point therefore eliminating much of the noise in the output signal [9,10,1 1], 


Noise generated by electrons 
not coming from photocathode 



Figure 2.6 Pulse height distribution of a) signal plus noise on a PMT 
b) PMT noise only. 

Good performance of a photon counting system requires that the SER pulses from 
the PMT remain narrow and that the discriminator is fast enough to respond to each 
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incoming input pulse. The PMT is load matched at 50 ohms with the amplifier to retain 
the narrow pulse width of the PMT output. The discriminator's ability to respond and 
recover to be ready for the next pulse is its pulse-pair resolution. A discriminator should 
have a short pulse-pair resolution in order to reduce pulse pile-up errors. Pulse pile-up 
happen at high light intensities when photon events occur near the same instance. If the 
discriminator cannot respond fast enough, these photon events may be counted as one 
pulse therefore causing an underestimation in the total counts. Pulse pile-up error e p can 
be estimated by Equation 2.8 [ 9 ], 

s p = l-exp (- 7 ■/?•/,) ( 2 . 8 ) 

where 7 is the QE of the PMT, R is the count rate of incoming pulses, and td is the 
discriminator's output pulse width. For a PMT with 7 = 31.5% and a discriminator pulse 
width of 5 nsec, the pulse pile-up error is less than 5% for a count rate of 30 MHz. 

2.3.3 Noise in photon counting systems 

Several types of noise exist in photon counting systems. These sources of noise 
come from the non-ideal characteristics of the photomultiplier tubes, background light, 
and standard electronic noise associated with the accompanying electronics. Each noise 
component reduces the signal-to-noise ratio, which limits the overall sensitivity of the 
detector 

The main source of noise in the PMT is called dark current. Dark current is the 
amount of current that results from electrons that are thermionically emitted by the 
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dynodes and photocathode when no incident light is present. Since the dark current is 
heavily temperature dependent, maintaining the PMT at a low temperature with liquid or 
air-cooling can reduce it. Other non-ideal mechanisms occur such as leakage current, 
glass scintillation, and field emission but are extremely small and can be avoided through 
proper PMT tube operation and handling precautions. It is usually assumed that the 
signal output of a PMT is linearly proportional to the input light intensity and when no 
light is present there is only the dark current present. While this is a good assumption 
for low light levels, the PMT output is nonlinear for high input light intensities. Also, 
when the PMT is exposed to momentary high light intensities, the output does not return 
to the dark current level immediately but instead decays slowly. This effect is called 
signal-induced noise (SIN), which can limit the range and accuracy of ozone 
measurements in a DIAL system [12,13,14,15,16], Characterization of signal-induced 
noise effects in DIAL systems can be found in Appendix A. 

Depending on the type of application, background light may be a large source of 
noise in photon counting. In lidar applications, background light from the atmosphere 
such as sunlight in the daytime or moonlight at night adds a DC baseline to the measured 
lidar signal. Narrowband optical filters are often used to reduce the background intensity. 
The remaining background counts can be measured separately when no lidar signal is 
being pulsed in the atmosphere. In the data analysis, these counts would then be 
subtracted from the measured lidar return signal counts. 

The final noise component in photon counting originates from the electronic 
components that make up the system. This type of noise is unavoidable but can be 
reduced in some situations. One of the best ways to reduce electronic noise is bv using 
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shielded coaxial cables that limit electromagnetic interference. Also, using cables that 
are short in length helps reduce electromagnetic interference, pulse broadening, and 
signal reflections. 

The number of pulses that can be counted in a given time can be described by 
Poisson statistics [9,10], If N ph is the average number of signal pulses, then the standard 
deviation known as shot noise is: 


n ph = (2.9) 

The total number of pulses received at the anode is a combination of pulses from the 
incident light signal, background light, and dark current. Since Nph cannot be measured 
separately, it must be calculated by measuring the background and dark current when no 
signal is present on the PMT and subtracting it from the total signal. The shot noise rib 
and rid from the background light and the dark current are 

(2.10a) 

n j = (2.10b) 

where N b and N d are the number of pulses from background and the dark current. 
Therefore, the total noise in the system nt ot can be calculated from each individual noise 
component as follows: 
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and substituting N ph = n ph 2 , N h = n b 2 , and Nj = nf 


(2.11) 
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The SNR is the ratio of the pulses from the signal to the total noise in the system. 
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The SNR in Equation 2. 13 is a function of the total number of counts measured in a given 
time [9,10], To calculate the SNR as a function of count rate (counts per second) for a 
measurement time t, we substitute N’ ph = {N ph /t), N' b = (N b /t ), and N’ d = (N d /t) 
and the SNR becomes 


SNR = 


•K»+2(lv; + Ar) 


(2.14) 


This is the famous result for photon counting. Equation 2.14 shows that the SNR for 
photon counting increases as the square root of the measurement time [9,10], Long 
counting times can reveal extremely low signals that would normally be lost using analog 
detection techniques. 

Photon counting offers a very sensitive method to detect low light levels 
that do not produce analog signals. The discriminator allows the elimination of electronic 
baseline effects that can cause errors in DIAL calculations. Also, the improvement in the 
SNR can increase the overall measurement range of a lidar system. 
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2 . 3.4 Photon Counting Form of the Lidar Equation 

The receiver system for this research effort uses photon counting as the light 
detection technique. The elastic backscattered lidar equation (Eq. 2.1) gives the power 
collected by a lidar receiver system from a laser pulsed into the atmosphere. This 
equation can be easily transformed from the average power received to the number of 
photons per second collected. The average power P a in the laser pulse is described as 



( 2 . 15 ) 


where E 0 is the energy in the laser pulse and r L is the laser pulse width. The number of 
photons n k contained in a laser pulse with energy E 0 is 


n, 


e p {X) 

A 


( 2 . 16 ) 


where E P (A ) is the energy of a photon at wavelength A, h is Planck's constant, and c is 
the speed of light. Substituting into Equation (2.1) yields the count rate C r (A,R) of 
photons collected by a lidar receiver system. 
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The photon counting form of the lidar equation (Eq. 2.8) will be used in Chapter 5 
'Results" to calculate predicted lidar signal returns. The backscatter and extinction terms 
will be obtained from Elterman 25 who has tabulated standard atmospheric attenuation 
coefficients for lidar applications. The predicted signals will be compared to actual lidar 
return signals measured with the fiber-optic coupled lidar receiver system. The DIAL 
equation, neglecting the backscatter and extinction terms, will be used to calculate ozone 
number densities as a function of altitude from the lidar return signals. 



CHAPTER 3 


DESIGN AND DEVELOPMENT OF THE 
FIBER-OPTIC COUPLED LIDAR RECEIVER 

This chapter describes the design of a fiber-optic coupled lidar receiver that will 
be used to measure the ozone concentrations in the stratosphere. The general design of 
the receiver system is discussed. The influence of each component is then covered to 
evaluate the performance characteristics. Finally, the total throughput efficiency of the 
receiver is calculated. 

3.1 FIBER-OPTIC COUPLED LIDAR RECEIVER 

The type of telescope configuration used in a lidar system depends on the 
application and cost considerations. Three popular types of telescopes, the Newtonian, 
Cassegranian, and the Gregorian, are shown in Figure 3. 1 [6], 



(ci Cassegramian 


Figure 3.1. Different receiver telescope configurations commonly used in lidar 
receivers. 
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The Newtonian configuration was chosen for this research because it offers the best way 
to couple light into the fiber-optic cable with little loss by obstruction of the primary 
mirror. Also, it is inexpensive and easy to construct. 

The receiver telescope for this research is intended for use with a UV DIAL 
system that transmits online and offline laser wavelengths of 300 and 311 run 
respectively. The purpose of the telescope is to collect as much of the backscattered lidar 
signal while neglecting as much of the background light as possible. The receiver system 
is divided into two parts, the telescope with fiber-optic output and the detector package as 
shown in Figure 3.2. 



Figure 3.2. Fiber-optic coupled receiver system for stratospheric ozone 

measurements 


The receiver telescope is made up of three optical components: the primary 
mirror, a small 45-degree turning mirror, and an optical fiber. The components, listed in 
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Table 3.1, are housed in a carbon-epoxy cylinder that is extremely strong and lightweight. 
The carbon-epoxy has low temperature sensitivity, and therefore does not expand and 
contract during environmental temperature changes. This helps avoid misalignment of 
the optics that can occur from material expansion and contraction. 

Light reflected from the primary mirror is focused onto the flat turning mirror 
and then focused onto the core of the optical fiber. With proper design all the light 
collected by the primary mirror can be coupled into the 1-mm diameter fiber-optic cable. 


Component 

Manufacturer 

Properties 

Primary 

Desmaris 

Parabolic Mirror 

Mirror 

Optical Systems 

355.6 mm diameter 


Inc. 

762 mm focal length 
aluminum with MgF3 coating 

Turning 

Newport 

25.4 mm dia. Flat mirror 

Mirror 

10D20.RM2 

Rmax = 98% @ 300 nm 45 
incidence 

Optical 

CeramOptec 

0.5 meter UV quartz 

Fiber 

UV 1000/1060 

numerical aperture = 0.28 
Core diameter = 1 mm 


Table 3.1 Optical components of fiber-optic coupled lidar receiver 

The purpose of the detector package is to filter and guide the light to the PMT 
detector. The output of the fiber is connected to the detector package. Light exits the 
fiber and is collimated by a 25.4-mm diameter plano-convex lens. The collimated light 
passes through the UV interference filter and strikes the 50.4 mm diameter photocathode 
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of the PMT. The electrical pulses from the PMT are connected to a photon counting 
system for data acquisition and analysis. Table 3.2 lists the components used in the 
detector package. 


Component 

Manufacturer 

Properties 

Collimating 

Lens 

Newport 

#SPX016.R2 

25.4 mm dia. 
Plano-Convex 
Focal length = 40 mm 

Interference Filter 

BARR# 

2290 

50% transmission at 300 nm 
40% transmission at 3 1 1 nm 
24 nm FWHM bandwidth 

PMT Detector 

Electron Tubes Inc. 
Type 9214Q#5150 

12 dynode chain linear focused 
50.4 mm dia. photocathode 
QE = 31.5% and 
Gain @1200 VDC =lxl0 7 


Table 3.2 Components of detector package. 


3.2 DESIGN CHARACTERISTICS OF THE RECEIVER 

The design characteristics of the receiver telescope require evaluation of its total 
throughput efficiency and its field of view (FOV). The total throughput efficiency 
depends upon the size, shape, alignment, and spectral reflectance/transmittance of each 
optical component. A high efficiency is desired in order to have a long ozone 
measurement range. The telescope FOV determines how much of the backscattered lidar 
signal is collected. The FOV of the telescope should be large enough to collect the 
backscattered laser signal but small enough to neglect unwanted atmospheric background 
light [17], 


T"7 









3.2.1 Primary Mirror 

The primary mirror is the main optical element in the telescope receiver. The 
diameter, reflectance, and focal spotsize determine how much of the lidar signal is 
collected and coupled into the fiber-optic cable. The amount of signal collected increases 
as the square of the diameter for the mirror. The primary mirror for the fiber-optic 
coupled receiver has a diameter of 35.5 cm and is coated with a standard MgFj plus 
aluminum coating for high reflection (90%) at UV wavelengths. The fiber-optic cable 
and the turning mirror are mounted on a 1 .225 cm wide aluminum bar that transverses the 
diameter of the telescope as shown in Figure 3.3. 



Figure 3.3 Front view of telescope showing fiber-optic cable and turning mirror 
mount blocking a portion of the primary mirror. 
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The total area of the mount blocking the primary mirror is 35.5 cm 2 . The area of the 
primary mirror is 993.1 cm 2 , therefore the total amount of light that will reach the 
primary is 


T 

mount 


993.1 - 35 .5 
993.1 


x 100 = 96.4% 


(3.1) 


Light reflected from the primary mirror is focused to a finite spotsize. The 
minimum diameter of the spotsize of a mirror is its diffraction-limited spotsize [18]. The 
focal point of a uniformly illuminated mirror has a minimum spotsize diameter d s 
described by 



d M 


(3.2) 


where fy and d\( are the focal length and diameter of the mirror respectively [18]. 
Therefore, at 300 nm, the primary mirror from Table 3.1 will have a diffraction limited 
spotsize diameter of 1.28 um, which is approximately 800 times less than the 1 mm core 
diameter of the fiber-optic cable. 

3.2.2 Fiber-Optic Cable 

Light from the primary mirror is turned by the 45-degree mirror and focused onto 
the fiber-optic cable. The main goal in the design and construction of this receiver 
system was to obtain the highest efficiency of light transmission and signal-to-noise ratio 
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as possible. To achieve high efficiency, the receiver should contain as few optical 
components as possible to reduce losses to reflection, absorption, and scattering. This 
also reduces difficulties in alignment that is associated with using multiple optic mirrors 
and lens. The fiber-optic cable provides an easy means to guide the collected light from 
the telescope to the detector that is usually accomplished by using multiple optics. This 
section gives a brief description of optical fibers as waveguides and the loss mechanisms 
that occur. 

An optical fiber is a cylindrical waveguide that is made up of low loss glass 
materials. The central part of the fiber is known as the core and the region surrounding 
the core is called the cladding. The core and cladding provide the lightwave guiding 
feature of the fiber. If light enters the fiber at the proper acceptance angle, determined by 
the indices of refraction, then the light can propagate in the core through the length of the 
fiber by total reflection at the core-cladding boundary without losing energy to refraction 
into the cladding. 

For a wave to propagate into a waveguide it must undergo total internal reflection 
at the waveguide core-cladding boundary [19], A wave in a medium with refractive 
index of nj, incident on a boundary to another medium with refractive index ri 2 , is 
partially reflected in medium 1 and partially transmitted in medium 2. The angle of 
transmitted light into medium 2 can be found from Snell's law. 

n i sin(^) = n 2 sin(# 2 ) (3.3) 
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If ni is less than nz, and the angle of incidence of the wave is greater than some critical 
angle, then the wave will be totally reflected into medium 1, the core [19,20], The 
critical angle is the maximum angle that light in medium 1 can intercept medium 2 and be 
totally reflected is given by 


0 C = sin" 1 



(3-4) 


Once the condition for total internal reflection is met, light can propagate in the 
fiber with very little loss. The lightwaves travel through the fiber with different 
transverse spatial distributions known as modes [20], The core and cladding 
characteristics determine the number of modes that can be supported by the fiber. It can 
be shown that higher order modes have larger divergence angles and travel longer 
distances through the fiber [19], The lower order modes travel closer to the propagation 
axis and therefore travel shorter distances. This effect is known as modal dispersion. 
Modal dispersion consequently causes a spreading of the pulse width. This effect can be 
ignored for short length fibers as in the case for the receiver for this research [19], 

To efficiently couple light into a fiber, the incident light must enter within the 
acceptance angle of the fiber. The acceptance angle is defined as the maximum angle in 
which light can enter the core of the fiber and still propagate by internal reflection 
through the fiber. A lightwave in air (refractive index n a = 1), incident on a fiber at angle 
0 refracts into the core. If the angle of the refracted light into the core is less than the 
critical angle when it strikes the core-cladding boundary, then the light will undergo total 
internal reflection (see Figure 3.4). The acceptance angle can be found by deriving the 
angles of refraction governed by Snell's law [19], A constant known as the numerical 
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aperture (NA) relates the indices of refraction of the core and cladding of the fiber. The 
maximum acceptance angle a aee of a fiber can then be calculated as 

a a cc ~ sin"’(AC4) (3.5) 

where 


NA = ^(n\ -nl) (3.6) 


n 0 



Figure 3.4 Acceptance angle of an optical fiber. 


The fiber from Table 3.1 with a NA - 0.28 has an acceptance angle of 16.26°, which is 
greater than the 13.8° mirror illumination angle 0in C as shown in Figure 3.5. 



= tan 


-i 


t j 


= tan 


./ 177.8' 


762 


= 138° 


Figure 3.5 Angle of incident light on fiber. 
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Therefore, all the light from the mirror should be coupled into the fiber except for a 
small percentage that is reflected at the air-fiber boundary (4%)[20], This reflection loss 
is calculated from standard electromagnetic theory using Snell's law [20], If the fiber end 
is not antireflection coated, then the reflection coefficient at the boundary 
is^? « 4% [21,22] yielding a total coupling efficiency into the fiber of 96%. 

Once light is launched into the fiber, the light is guided to the output subject to 
losses in the fiber. There are three different types of loss mechanisms in optical fibers; 
the first is from bending of the fiber, the second is from scattering losses, and third are 
losses due to absorption [19,23], It is fairly obvious why bends in fibers should give rise 
to losses. To maintain a wavefront that is perpendicular to the direction of propagation 
when the fiber is bent, the part of the wave that is close to the core-cladding boundary 
must travel faster than the part near the center of the core. In theory, each mode extends 
infinitely into the cladding even though the electric field decreases exponentially as it 
projects into the cladding. That would mean when there is a bend in the fiber, the 
extended light on the outside of the bend would have to travel faster than the light at the 
inside of the bend in order to maintain a uniform wavefront. If that were the case, the 
outer part of the wavefront would have to travel faster than the speed of light to keep up, 
which according to relativistic theory is not possible. Therefore, the light not contained 
in the core will be lost to radiation (heat) from the core. Fibers with bends of small radii 
of curvature will cause higher order modes that are propagating in the fiber to incur the 
greatest loss [19,23], 
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The main losses resulting from the physical structure of the fiber are scattering 
losses and absorption losses. The scattering losses occur from the chemical structure of 
the fiber and the lattice shape of the molecules. Ideally the lattice structure of the 
molecules would be perfectly placed and symmetric throughout the material, but in 
reality this is not the case. The molecules are connected in a rather random fashion [23], 
The chemical composition can also vary from place to place throughout the fiber. Both 
of these types of non-ideal structure and chemical composition can cause differences in 
the refractive indices of the material. If these differences are on the order of the 
wavelength divided by ten or more, then this fluctuation can act as Rayleigh scattering 
points. The Rayleigh scattering coefficient is proportional to MX* and is the cause of 
most of the attenuation at shorter wavelengths [23], To reduce the loss at UV 
wavelengths by Rayleigh scattering, materials such as UV grade fused silica must be 
used. 

Losses in the longer wavelengths result from two absorption processes in the fiber 
material. The first process is a result of the transitions between lattice vibrational states. 
These transitions correspond to the infrared region from about 1600nm and above. 
Impurities in the fiber material are the other source of attenuation. Often, the wavelength 
of the transmitted light corresponds to a peak absorption wavelength of the impurity. 
One of the largest contributors of this type of loss is the presence of hydroxyl ions (OH-). 
The major advances in decreasing the attenuation in fiber optics over the years have come 
from the effective control over the densities of impurities during the manufacturing 
process [23], 
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The fiber used in this research effort (Table 3.1) is made of UV grade fused silica 
(quartz) that has a high transmission at UV wavelengths with -10 dB/meter loss at 300 
nm. The fiber is 0.5 m long giving a total loss of only 1%. The coupling efficiency, loss 
by reflection at each end, and the transmission of the fiber yield a high total throughput 
efficiency of the fiber to be 91%. Therefore, less than 10% of the total signal should be 
lost using the fiber-optic cable to guide the signal to the detector package. 

3.2.3 Detector Package 

Optimum performance of the entire receiver requires that the filter passes only the 
desired wavelengths and that all of this light strikes the photocathode of the PMT. 
Interference filters are best suited for UV wavelengths because of their high UV 
transmission. Since the spectral transmission of interference filters is dependent on the 
incident angle, the light from the output of the fiber must be collimated before being 
passed through the filter. 

The collimating lens has a diameter of 25.4 mm and is anti-reflection coated to 
provide > 99% transmission at 300 nm. The fiber is placed at the focal length of the lens 
to collimate the light exiting the fiber. The spectral transmission of the interference filter 
was measured at the NASA Langley spectroscopy lab with a Cary-14 spectrophotometer. 
The UV interference filter has a total bandwidth of 24 nm FWHM from 288 to 312 nm. 
The transmission at the laser wavelengths 300 and 3 1 1 nm is 50% and 40% respectively. 

The collimated light passes through the filter and strikes the photocathode of the 
PMT. ). The diameter of the photocathode is 50.4 mm, therefore all of the collimated 
light that is passed through the filter strikes the photocathode. The PMT used in this 
research effort was an £77 9214Q. This PMT was suggested for DIAL applications by 
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Eccles 26 based on its high QE at 300 nm (31.5%) and its high gain (1 x 10 7 ) at low 
operating voltage of -1200 VDC (see Appendix c) 

3 . 2.4 Telescope FOV and Total Throughput Efficiency 

The FOV of the telescope is limited by the size of the field stop aperture. In the 
case of the fiber-optic coupled receiver, the core of the fiber is the field stop. At short 
ranges, the FOV of each point of the mirror is slightly different [6], This effect decreases 
for longer ranges and can be ignored for this receiver. The application of this receiver 
system is to measure upper atmospheric (stratosphere) ozone. At these high altitudes, the 
full angle FOV can be estimated as 

d f 

®FOV ~ r (3 6 ) 

J u 

where df is the field stop diameter and/vr is the focal length of the primary mirror [6], 
The FOV of a telescope should be large enough to collect all the backscattered radiation 
while minimizing the amount of background light collected. For maximum signal-to- 
noise ratio, the FOV of a receiver (depending on the application) should be between 1.0 
and 2.0 times the divergence of the laser beam [17]. The lasers for NASA Langley's UV 
DIAL system typically have divergences between 0.5 and 1 mrad [7], The FOV of the 
fiber-optic coupled lidar receiver, with the 762-mm focal length primary mirror and 1- 
mm diameter fiber core, is 1.3 mrad. This FOV is within the constraints defined in the 
Lidar Technical Notes 17 . 

The total throughput efficiency of the receiver telescope to the PMT photocathode 
can be calculated from the individual spectral characteristics of each part of the system. 
The total efficiency is then 
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where Tm,,,,,,, is the transmission of the light with the cable and turning mirror mount 
blocking the primary, Rmi(A) is the reflectivity of the primary mirror, R\a( A ) is the 
reflectivity of the turning mirror Cfiber is the coupling efficiency of light into the fiber, 
T f {X) is the transmission of the fiber, Tc 0 u( A) is the transmission of the collimating lens, 

and T[f{ A) is the transmission of the interference filter at wavelength A . Therefore, the 
expected total throughput efficiency for the fiber-optic coupled lidar receiver to the PMT 
detector is 

s( A) =)0.964 • 0.90 • 0.98 • 0.96 • 0.95-0.99 • 0.50 = 0.387 = 38.7% @ 300 nm 
s{A) = 0.964 • 0.90 • 0.98 • 0.96 • 0.95 • 0.99 • 0.40 = 0.307 = 30.7% @ 3 1 1 nm 

The total efficiency of the telescope receiver calculated above does not include the QE 
efficiency of the PMT detector. We can solve for the system transmission coefficient 
^{A) in the photon counting form of the lidar equation by multiplying QE of the detector 
by the total throughput efficiency. If we assume that all electron pulses produced by 
incident photons are counted, then the system transmission coefficient g(A) is 

c(A) = s{A)-OE (3.8) 
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The performance characteristics of the entire receiver system are expected to yield 
a high signal-to-noise ratio. This system has a high transmission at the laser wavelengths 
of 300 and 3 1 1 nm, and a low transmission at wavelengths above and below this region. 
The above characteristics should allow all the backscattered laser energy, where there is 
full overlap of the FOV and laser beam, to be collected while neglecting much of the 
background light that may be present. 
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CHAPTER 4 


EXPERIMENTAL SETUP 

The following chapter describes the experimental setups for measuring the 
spectral transmission characteristics of the telescope and detector package. The current 
ozone DIAL system at NASA Langley Research center is also described. Finally, the 
experimental setup for using the fiber-optic coupled receiver system to make ground 
based stratospheric ozone measurements is discussed. 

4.1 SPECTRAL TRANSMISSION CHARACTERISTICS OF RECEIVER 
SYSTEM 

The performance of the fiber-optic coupled telescope depends upon the optics that 
make up the receiver. The spectral characteristics of each optical component must be 
measured in order to determine the expected signal-to-noise ratio of the system. High 
signal-to-noise ratio translates into long atmospheric measurement ranges, which is 
desired for remote sensing of ozone in the high altitude stratosphere. 

In the experimental setup, a collimated light source was used as shown in Figure 
4.1. A 1 000- W Xe continuous wave lamp placed in an air-cooled housing provided an 
intense UV light spectrum. The lamp was powered by an adjustable DC power supply 
that has an optical feedback amplifier coupled into the lamp housing. The feedback 
amplifier provides lamp stability to under 0. 1% ripple. A high intensity grating 
monochrometer was used to select the desired wavelengths from the Xe lamp. Light 
from the monochrometer was coupled into a UV enhanced quartz fiber with 1 mm core 
diameter. The output of the fiber acted as a point source that was placed at the focal 
point of a collimating mirror. An adjustable iris was placed at the output of the fiber to 





control the diameter of the illumination circle on the collimating mirror. The cone of 
light exiting the fiber was collimated into a parallel beam that was reflected into the fiber- 
optic telescope designed in this research effort. 


Collimating Mirror 



Figiire 4.1. Collimated light source for measuring spectral characteristics of 

receiver optical elements. 


The spectral response of the telescope primary mirror, turning mirror, and fiber 
were measured by placing a standard UV PIN diode detector at each point as shown in 
Figure 4.2. The current from the detector was measured with a picoammeter. The 







reflectivity or transmission of each element was determined by taking the ratio of the 
measured reflected/transmitted light and the measured incident light. 
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Figure 4.2 Measurement points A,B,C,D and E for spectra] transmission and 
reflection of each optical element. 

The collimated light source strikes the primary mirror, which focuses the light to a 
finite spotsize at its focal length. The PIN diode detector was first placed at the output of 
the light source fiber (see Figure 4.1) to measure the incident light on the telescope. This 
measurement was multiplied by the reflectivity of the collimating mirror in order to 
calculate the incident light striking the primary mirror. The detector was then placed at 
point A near the turning mirror. In order to make as accurate measurement as possible, 
the full area of the detector was used. The detector was then placed at point B to measure 



the reflectivity of the turning mirror. The measurements from point A were used as the 
incident light on the turning mirror. 

The coupling efficiency and the absorption of the fiber were measured by taking 
the ratio of the output and input light intensities at points B and C respectively. This ratio 
yielded the coupling efficiency of the fiber that includes the amount of light coupled into 
the fiber with losses by reflection at the endpoints and absorption and scattering losses as 
the light travels through the fiber. 

The previous measurement did not include losses due to the bending of the fiber. 
The fiber in the configuration for this receiver telescope has one 90-degree bend with a 
radius of curvature of 17. 1 cm. The measurement of the output of the fiber with the bend 
at point D divided by that of the measurement without the bend at point C is the amount 
of bending loss. Losses from multiple bends in the fiber were not measured because the 
configuration of the telescope requires only one bend in the fiber before the output is 
placed into the detector package. 

Since the spectral characteristics of the interference filter and the PMT are known 
(see Ch j. 2.3), the collimating lens is the only component in the detector package to be 
measured. The transmission of the collimating lens was measured by taking the ratio of 
the incident light measured at point D and the transmitted light at point E. The area of 
the PIN diode detector is only 3 cm 2 so the fiber was placed a few centimeters off of the 
focal point. This created an illumination circle that could be measured by the full area of 


the PIN diode detector. 



4.2 STRATOSPHERIC OZONE MEASUREMENT SYSTEM SETUP 
4.2. 1 Ozone DIAL Laser System 

The current ozone DIAL laser system at NASA Langley Research Center, shown 
in Figure 4.3, consists of two types of lasers to produce the online and offline 
wavelengths. Two high-power Continuum Nd:YAG lasers are frequency-doubled, which 
then pump dye lasers. The dye lasers are doubled to produce 289 and 299 nm for 
tropospheric measurements or 301 and 311 nm for stratospheric measurements. These 
lasers operate at a repetition rate of 30 Hz, The output power depending on the type of 
dye used is between 10 and 30 mJ per pulse. The laser beams are transmitted in the 
zenith for ground based measurements through a door in the building roof {Richter], 

The DIAL system has two 35.5 cm diameter Cassegrainian telescopes to collect 
the returned light backscattered from the atmosphere from the zenith and nadir directions. 
Each optical surface is protected with a triple-band anti-reflection coating to reduce 
losses of the light passing through the receiver optics. 

Photomultiplier tubes (PMTs) are used in an analog configuration for light 
detection. A 5 MHz analog-to-digital converter (ADC) connected to a computer stores an 
average of a variable number of lidar returns. The 5 MHz ADC allows a 30-meter 
atmospheric measurement resolution. 

The ozone DIAL system was used as the laser transmitter source for this research. 
The fiber-optic coupled lidar receiver using photon counting was used to measure the 
backscattered lidar returns from the transmitted online and offline laser pulses. No 
measurements were made by the NASA analog system because the beams were aligned 
to the fiber-optic coupled receiver telescope. 
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Figure 4.3. Current ozone UV DL\L system in ground test lab. Fiber-optic coupled 
receiver system placed one meter away from lasers (drawing not to 
scale). 

4.2.2 Fiber-optic coupled receiver system 

Stratospheric ozone measurements were made using the fiber-optic coupled 
receiver shown in Figure 4.4. The receiver system was connected to the photon counting 
system to measure the backscattered lidar signals from the online and offline lasers. Data 
from the photon counting system was sent to a computer for analysis. 
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The fiber-optic coupled lidar feceivw Was- placed approximately one meter away 
from the exit of the laser beams on an adjustable tripod (see Figure 4.3). The 
micrometers allowed for additional adjustments that were sometimes needed to align the 
telescope and laser beams. This biaxial arrangement allowed discrimination against 
strong near-field returns that saturate the PMT detector. The return signals were 
maximized by monitoring the output from the PMT anode on an oscilloscope. After the 
telescope and the lasers were aligned, the anode of the PMT was connected to the photon 
counting system to measure the backscattered lidar returns (see Fig. 4.5). 



The photon counting system used in this research effort is shown in Figure 4.5. 
Although the receiver telescope was placed one meter away from the lasers, some of the 
strong near-field signal was still in the FOV of the telescope. Therefore, the PMT is 
gated with the circuit in Appendix C to avoid saturation. 
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Figure 4.5 Photon counting system setup, 




Signal pulses from the PMT were jsent to a Standford Research Systems 300 MHz 
amplifier with a gain of 5. The amplifier has an input impedance of 50 ohms so the PMT 
is properly load matched. The amplified pulses were then passed through a 300 MHz 
Phillips Scientific 704 Quad discriminator. The discriminator has a variable threshold 
setting from -10 mV to -1 V and a variable output pulse width that ranges from 2 nsec to 
50 nsec. The input pulses above the set threshold were sent to a DSP 2090S 
multichannel-scaler (MCS), When the MCS receives a trigger signal, it counts the 
number of pulses received from the discriminator and stores the number into time bins. 
The counting time for each bin is set by the clock rate. A 1 MHz clock rate yields a 
count time of 1 usee per bin. The number of pulses in each bin is stored in a DSP 4 101 A 
averaging memory. The 1 MHz clock and 30 Hz laser trigger were provided by the 
NASA UV DIAL Transmitter Master Controller. 

The optimum operating conditions for these components in the photon counting system 
were determined in previous research by Eccles 26 . His research suggested that the PMT 
in Table 4. 1 be operated at a voltage of 1200 VDC with a four dynode gating circuit for 
maximum linear signal response. Eccles 26 also concluded that the discriminator output 
pulse width be set to 5 nsec and the discriminator threshold be set to about 110 mV [see 
Appendix CJ. 

The discriminator, MCS, and averaging memory are mounted in a standard 
CAMAC crate. A LeCroy 8901 A GPIB CAMAC crate controller allows computer 
interface with the mounted components. A computer with a National Instruments AT- 
GP IB/TNT interface allows control and data acquisition processes with the CAMAC 


crate. 



4.2.3 Computer Software for photon counting system 

A software program was developed in Microsoft™ VISUAL C++ to control and 
acquire data from the photon counting system. The program is a 32-bit WINDOWS™ 
application that provides ease of operation and a user-friendly environment to control the 
photon counting system. The program records setup information from the user including 
laser repetition rate, laser pulse width, discriminator threshold voltage, discriminator 
output pulse width, and PMT operating voltage. The program sets the number of bins 
and the total photon counting time to accumulate data. After the total counting time has 
expired, the data from the averaging memory is read into the computer through the GPEB 
communication interface. The setup information and the data are recorded into 
individual files and displayed on the screen monitor. The program also includes a DIAL 
calculation feature that determines real time ozone conc entrations as a -function of altitude 
from the data acquired. For post data analysis, the program also converts the photon 
counting data files into text files that can be read by standard spreadsheet and 
mathematical software applications. 
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